3.16 \(\int \frac {1}{a+b \cosh ^2(x)} \, dx\)

Optimal. Leaf size=29 \[ \frac {\tanh ^{-1}\left (\frac {\sqrt {a} \tanh (x)}{\sqrt {a+b}}\right )}{\sqrt {a} \sqrt {a+b}} \]

[Out]

arctanh(a^(1/2)*tanh(x)/(a+b)^(1/2))/a^(1/2)/(a+b)^(1/2)

________________________________________________________________________________________

Rubi [A]  time = 0.02, antiderivative size = 29, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, integrand size = 10, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.200, Rules used = {3181, 208} \[ \frac {\tanh ^{-1}\left (\frac {\sqrt {a} \tanh (x)}{\sqrt {a+b}}\right )}{\sqrt {a} \sqrt {a+b}} \]

Antiderivative was successfully verified.

[In]

Int[(a + b*Cosh[x]^2)^(-1),x]

[Out]

ArcTanh[(Sqrt[a]*Tanh[x])/Sqrt[a + b]]/(Sqrt[a]*Sqrt[a + b])

Rule 208

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-(a/b), 2]*ArcTanh[x/Rt[-(a/b), 2]])/a, x] /; FreeQ[{a,
b}, x] && NegQ[a/b]

Rule 3181

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]^2)^(-1), x_Symbol] :> With[{ff = FreeFactors[Tan[e + f*x], x]}, Dist
[ff/f, Subst[Int[1/(a + (a + b)*ff^2*x^2), x], x, Tan[e + f*x]/ff], x]] /; FreeQ[{a, b, e, f}, x]

Rubi steps

\begin {align*} \int \frac {1}{a+b \cosh ^2(x)} \, dx &=\operatorname {Subst}\left (\int \frac {1}{a-(a+b) x^2} \, dx,x,\coth (x)\right )\\ &=\frac {\tanh ^{-1}\left (\frac {\sqrt {a} \tanh (x)}{\sqrt {a+b}}\right )}{\sqrt {a} \sqrt {a+b}}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]  time = 0.07, size = 29, normalized size = 1.00 \[ \frac {\tanh ^{-1}\left (\frac {\sqrt {a} \tanh (x)}{\sqrt {a+b}}\right )}{\sqrt {a} \sqrt {a+b}} \]

Antiderivative was successfully verified.

[In]

Integrate[(a + b*Cosh[x]^2)^(-1),x]

[Out]

ArcTanh[(Sqrt[a]*Tanh[x])/Sqrt[a + b]]/(Sqrt[a]*Sqrt[a + b])

________________________________________________________________________________________

fricas [B]  time = 0.46, size = 293, normalized size = 10.10 \[ \left [\frac {\log \left (\frac {b^{2} \cosh \relax (x)^{4} + 4 \, b^{2} \cosh \relax (x) \sinh \relax (x)^{3} + b^{2} \sinh \relax (x)^{4} + 2 \, {\left (2 \, a b + b^{2}\right )} \cosh \relax (x)^{2} + 2 \, {\left (3 \, b^{2} \cosh \relax (x)^{2} + 2 \, a b + b^{2}\right )} \sinh \relax (x)^{2} + 8 \, a^{2} + 8 \, a b + b^{2} + 4 \, {\left (b^{2} \cosh \relax (x)^{3} + {\left (2 \, a b + b^{2}\right )} \cosh \relax (x)\right )} \sinh \relax (x) - 4 \, {\left (b \cosh \relax (x)^{2} + 2 \, b \cosh \relax (x) \sinh \relax (x) + b \sinh \relax (x)^{2} + 2 \, a + b\right )} \sqrt {a^{2} + a b}}{b \cosh \relax (x)^{4} + 4 \, b \cosh \relax (x) \sinh \relax (x)^{3} + b \sinh \relax (x)^{4} + 2 \, {\left (2 \, a + b\right )} \cosh \relax (x)^{2} + 2 \, {\left (3 \, b \cosh \relax (x)^{2} + 2 \, a + b\right )} \sinh \relax (x)^{2} + 4 \, {\left (b \cosh \relax (x)^{3} + {\left (2 \, a + b\right )} \cosh \relax (x)\right )} \sinh \relax (x) + b}\right )}{2 \, \sqrt {a^{2} + a b}}, \frac {\sqrt {-a^{2} - a b} \arctan \left (\frac {{\left (b \cosh \relax (x)^{2} + 2 \, b \cosh \relax (x) \sinh \relax (x) + b \sinh \relax (x)^{2} + 2 \, a + b\right )} \sqrt {-a^{2} - a b}}{2 \, {\left (a^{2} + a b\right )}}\right )}{a^{2} + a b}\right ] \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*cosh(x)^2),x, algorithm="fricas")

[Out]

[1/2*log((b^2*cosh(x)^4 + 4*b^2*cosh(x)*sinh(x)^3 + b^2*sinh(x)^4 + 2*(2*a*b + b^2)*cosh(x)^2 + 2*(3*b^2*cosh(
x)^2 + 2*a*b + b^2)*sinh(x)^2 + 8*a^2 + 8*a*b + b^2 + 4*(b^2*cosh(x)^3 + (2*a*b + b^2)*cosh(x))*sinh(x) - 4*(b
*cosh(x)^2 + 2*b*cosh(x)*sinh(x) + b*sinh(x)^2 + 2*a + b)*sqrt(a^2 + a*b))/(b*cosh(x)^4 + 4*b*cosh(x)*sinh(x)^
3 + b*sinh(x)^4 + 2*(2*a + b)*cosh(x)^2 + 2*(3*b*cosh(x)^2 + 2*a + b)*sinh(x)^2 + 4*(b*cosh(x)^3 + (2*a + b)*c
osh(x))*sinh(x) + b))/sqrt(a^2 + a*b), sqrt(-a^2 - a*b)*arctan(1/2*(b*cosh(x)^2 + 2*b*cosh(x)*sinh(x) + b*sinh
(x)^2 + 2*a + b)*sqrt(-a^2 - a*b)/(a^2 + a*b))/(a^2 + a*b)]

________________________________________________________________________________________

giac [A]  time = 0.13, size = 39, normalized size = 1.34 \[ \frac {\arctan \left (\frac {b e^{\left (2 \, x\right )} + 2 \, a + b}{2 \, \sqrt {-a^{2} - a b}}\right )}{\sqrt {-a^{2} - a b}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*cosh(x)^2),x, algorithm="giac")

[Out]

arctan(1/2*(b*e^(2*x) + 2*a + b)/sqrt(-a^2 - a*b))/sqrt(-a^2 - a*b)

________________________________________________________________________________________

maple [B]  time = 0.09, size = 78, normalized size = 2.69 \[ -\frac {\ln \left (\sqrt {a +b}\, \left (\tanh ^{2}\left (\frac {x}{2}\right )\right )-2 \sqrt {a}\, \tanh \left (\frac {x}{2}\right )+\sqrt {a +b}\right )}{2 \sqrt {a}\, \sqrt {a +b}}+\frac {\ln \left (\sqrt {a +b}\, \left (\tanh ^{2}\left (\frac {x}{2}\right )\right )+2 \sqrt {a}\, \tanh \left (\frac {x}{2}\right )+\sqrt {a +b}\right )}{2 \sqrt {a}\, \sqrt {a +b}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(a+b*cosh(x)^2),x)

[Out]

-1/2/a^(1/2)/(a+b)^(1/2)*ln((a+b)^(1/2)*tanh(1/2*x)^2-2*a^(1/2)*tanh(1/2*x)+(a+b)^(1/2))+1/2/a^(1/2)/(a+b)^(1/
2)*ln((a+b)^(1/2)*tanh(1/2*x)^2+2*a^(1/2)*tanh(1/2*x)+(a+b)^(1/2))

________________________________________________________________________________________

maxima [B]  time = 0.42, size = 53, normalized size = 1.83 \[ -\frac {\log \left (\frac {b e^{\left (-2 \, x\right )} + 2 \, a + b - 2 \, \sqrt {{\left (a + b\right )} a}}{b e^{\left (-2 \, x\right )} + 2 \, a + b + 2 \, \sqrt {{\left (a + b\right )} a}}\right )}{2 \, \sqrt {{\left (a + b\right )} a}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*cosh(x)^2),x, algorithm="maxima")

[Out]

-1/2*log((b*e^(-2*x) + 2*a + b - 2*sqrt((a + b)*a))/(b*e^(-2*x) + 2*a + b + 2*sqrt((a + b)*a)))/sqrt((a + b)*a
)

________________________________________________________________________________________

mupad [B]  time = 0.35, size = 267, normalized size = 9.21 \[ -\frac {\mathrm {atan}\left (\frac {b^2\,{\mathrm {e}}^{2\,x}\,{\left (-a^2-b\,a\right )}^{3/2}\,\left (\frac {4\,\left (4\,a+2\,b\right )\,\left (8\,a^3+12\,a^2\,b+4\,a\,b^2\right )}{b^5\,{\left (-a^2-b\,a\right )}^{3/2}\,\sqrt {-a\,\left (a+b\right )}}+\frac {2\,\left (8\,a^2+8\,a\,b+b^2\right )\,\left (8\,a^2\,\sqrt {-a^2-b\,a}+b^2\,\sqrt {-a^2-b\,a}+8\,a\,b\,\sqrt {-a^2-b\,a}\right )}{a\,b^5\,\left (a+b\right )\,{\left (-a^2-b\,a\right )}^{3/2}}\right )}{4}+\frac {\left (2\,a^2\,b+2\,a\,b^2\right )\,\left (4\,a+2\,b\right )}{b^3\,\sqrt {-a\,\left (a+b\right )}}+\frac {\left (b^2\,\sqrt {-a^2-b\,a}+2\,a\,b\,\sqrt {-a^2-b\,a}\right )\,\left (8\,a^2+8\,a\,b+b^2\right )}{2\,a\,b^3\,\left (a+b\right )}\right )}{\sqrt {-a^2-b\,a}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(a + b*cosh(x)^2),x)

[Out]

-atan((b^2*exp(2*x)*(- a*b - a^2)^(3/2)*((4*(4*a + 2*b)*(4*a*b^2 + 12*a^2*b + 8*a^3))/(b^5*(- a*b - a^2)^(3/2)
*(-a*(a + b))^(1/2)) + (2*(8*a*b + 8*a^2 + b^2)*(8*a^2*(- a*b - a^2)^(1/2) + b^2*(- a*b - a^2)^(1/2) + 8*a*b*(
- a*b - a^2)^(1/2)))/(a*b^5*(a + b)*(- a*b - a^2)^(3/2))))/4 + ((2*a*b^2 + 2*a^2*b)*(4*a + 2*b))/(b^3*(-a*(a +
b))^(1/2)) + ((b^2*(- a*b - a^2)^(1/2) + 2*a*b*(- a*b - a^2)^(1/2))*(8*a*b + 8*a^2 + b^2))/(2*a*b^3*(a + b)))
/(- a*b - a^2)^(1/2)

________________________________________________________________________________________

sympy [A]  time = 46.72, size = 12026, normalized size = 414.69 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(a+b*cosh(x)**2),x)

[Out]

Piecewise((zoo*tanh(x/2)/(tanh(x/2)**2 + 1), Eq(a, 0) & Eq(b, 0)), (2*tanh(x/2)/(b*(tanh(x/2)**2 + 1)), Eq(a,
0)), (-tanh(x/2)/(2*b) - 1/(2*b*tanh(x/2)), Eq(a, -b)), (-5*I*a**(5/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a +
b) + a/(a + b) - b/(a + b))*log(-sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a
**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) +
a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt
(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) -
b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a
+ b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-
2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
+ 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b))) + 5*I*a**(5/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(sqrt
(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqr
t(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)
*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a +
b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/
(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2
*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a +
b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)
/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - 3*I*a**(5/2)*sq
rt(b)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(-sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)
)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/
(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*
sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 1
0*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a +
b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqr
t(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sq
rt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) + 3*I*a**(5/2)*sqrt(b)*sqrt(2*I*sqrt(a)*sqrt(b)/(a +
b) + a/(a + b) - b/(a + b))*log(sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a
**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) +
a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt
(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) -
b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a
+ b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-
2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
+ 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b))) + 10*I*a**(3/2)*b**(3/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(-s
qrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*
sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3
/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(
a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(
b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqr
t(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt
(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - 10*I*a**(3/2
)*b**(3/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a
/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt
(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-
2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
- 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(
a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)
*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)
)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - 2*I*a**(3/2)*b**(3/2)*sqrt(2*I*sqrt(a)*sqrt(b)/
(a + b) + a/(a + b) - b/(a + b))*log(-sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/
(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a +
b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)
)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt
(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*
sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b)
+ a/(a + b) - b/(a + b))) + 2*I*a**(3/2)*b**(3/2)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*lo
g(sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt
(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a
**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) +
a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*s
qrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) +
a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*
sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - I*sqrt(a
)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(-sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) +
a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqr
t(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(
-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)
) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/
(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a
)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b
))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) + I*sqrt(a)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(
a + b) + a/(a + b) - b/(a + b))*log(sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*
I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b)
+ a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*s
qrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b
) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)
/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqr
t(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a +
b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a
/(a + b) - b/(a + b))) + I*sqrt(a)*b**(5/2)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(-sqr
t(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*s
qrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/
2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b
)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt
(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(
b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - I*sqrt(a)*b**
(5/2)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a +
b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(
a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*s
qrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10
*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b
) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt
(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqr
t(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - a**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) -
b/(a + b))*log(-sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*
sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqr
t(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt
(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b
) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt
(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt
(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b
))) + a**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a
/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt
(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-
2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
- 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(
a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)
*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)
)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - a**3*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a +
b) - b/(a + b))*log(-sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqr
t(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) -
b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a
)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(
a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)
*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3
*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(
a + b))) + a**3*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b
) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) -
b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)
*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*s
qrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b)
+ a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sq
rt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) + 10*a**2*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b
) + a/(a + b) - b/(a + b))*log(-sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a*
*(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a
/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(
2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) -
b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a
+ b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2
*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
+ 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b))) - 10*a**2*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(sqrt(2*I*sqrt(a)
*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b)
+ a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sq
rt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a +
b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/
(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*s
qrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a +
b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a
/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) + 2*a**2*b*sqrt(2*I*sqrt(a)*s
qrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(-sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh
(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(
b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/
(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b)
+ a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(
a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**
2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b)
- b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(
a + b) + a/(a + b) - b/(a + b))) - 2*a**2*b*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(sqrt
(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sq
rt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2
)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)
/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(
2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a +
b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b
)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - 5*a*b**2*sqrt(
-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(-sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*
sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a
/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqr
t(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sq
rt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a +
b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b
) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt
(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) + 5*a*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a
+ b))*log(sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2
*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))
- 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a
+ b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sq
rt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(
a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a
+ b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*s
qrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) + 3
*a*b**2*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*log(-sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(
a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sqrt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a +
b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b
)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*
I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) -
10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*s
qrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*
sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))) - 3*a*b**2*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a
+ b) - b/(a + b))*log(sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + tanh(x/2))/(8*I*a**(7/2)*sq
rt(b)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b)
- b/(a + b)) - 8*I*a**(3/2)*b**(5/2)*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(
a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a**4*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)
)*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**3*b*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/
(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) - 10*a**2*b**2*sqrt(-2*I*sqrt(a
)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b)) + 2*a*b**
3*sqrt(-2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/(a + b))*sqrt(2*I*sqrt(a)*sqrt(b)/(a + b) + a/(a + b) - b/
(a + b))), True))

________________________________________________________________________________________